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Motivated by recent experimental and numerical results, a simple unifying picture of intermittency 
in turbulent shear flows is suggested. Integral Structure Functions (ISF), taking into account explic- 
itly the shear intensity, are introduced on phenomenological grounds. ISF can exhibit a universal 
scaling behavior, independent of the shear intensity. This picture is in satisfactory agreement with 
both experimental and numerical data. Possible extension to convective turbulence and implication 
on closure conditions for Large-Eddy Simulation of non-homogeneous flows are briefly discussed. 



PACS: 47.27-i, 47.27.Nz, 47.27.Ak 

Statistical properties of turbulent flows are usually 
characterized in terms of the scaling behavior of velocity 
Structure Functions (SF). These quantities are defined 
as the statistical moments of longitudinal velocity incre- 
ments across a separation r at the location x: D p (x, r) — 
(Sv(x,r) p ). In homogeneous and isotropic turbulence, 
D p (x,r) only depends on the distance (or scale) r. Ex- 
perimental and numerical observations support the idea 
that D p (r) display universal power-law dependence on r 
in the so-called inertial range, i.e. D p {r) ~ r^ v . Uni- 
versality refers here to the scaling exponents £ p being 
independent of the stirring process of turbulence. The Q p 
values are found to be in disagreement with Kolmogorov's 
linear prediction (K41) ( p = p/3 fm. The understanding 
of this correction to K41, usually referred to as intermit- 
tency, has stimulated many phenomenological and theo- 
retical works during the last 40 years (see Q for a recent 
review). 

Only recently, interests in understanding intermit- 
tency in non-homogeneous turbulent flows have started 
to emerge (see p|-|Tc|). The major point is to understand 
how the phenomenology of intermittency is modified, or 
can be extended, in case of non-homogeneous flows. A 
common characteristic of such flows, e.g. wall-bounded 
flows, is the presence of a non-zero mean velocity gradient 
(usually called shear). Note that shear does not neces- 
sarily imply inhomogeneity, e.g. in homogeneous shear, 
straining and rotational flows. 

Our investigation starts from the following key obser- 
vations: i) In presence of a strong shear, intermittency, 
defined as the deviation of scaling exponents £ p from the 
linear law, is larger than in homogeneous and isotropic 
turbulence, ii) Relative scaling exponents, measured in 
very different flows but in positions where the shear is 
strong enough, seems to be very similar (universal). 

Data in Fig. [j] fully confirm our two key observations. 
They come from very different situations: near the wall in 
a channel flow numerical simulation and experiment 
p] , in the logarithmic sublayer of a boundary layer flow 



[J7j] , near a strong vortex [g| 
and in a Kolmogorov flow 



in the wake of a cylinder 
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FIG. 1. Intermittency corrections to scaling exponents, 
~ p/3, in homogeneous and isotropic turbulence (□), 
and several turbulent shear flows (A). Intermittency correc- 
tions are significantly larger in presence of shear and display 
universality. 

In order to provide a theoretical unifying framework 
for non-homogeneous turbulent shear flows, we start from 
the Navier-Stokes (N-S) equations. The velocity field can 
be decomposed into a mean value (average is meant on 
time) plus a fluctuating part: v(x;t) = v(x) + v'(x;t). 
It yields the usual Reynolds decomposition 



D f v' 



- S i3 {x)v' 3 + v'jdjv'i - v'jdjvl = -d^ + uAv'i (1) 



with D t = (d t + Vjdj). The shear is defined as Sij (x) — 
djVi (x) and will depend on the mean flow geometry. In 
regions where Sij — 0, e.g. very far from the boundaries, 
turbulence can be considered as homogeneous. Other- 
wise, the shear term &. 



(x)v'j must be taken into account. 



The results displayed in Fig. [I] show that the presence of 
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this term modifies significantly the statistical properties 
of turbulence. In order to better hilight the physical im- 
plication of the shear, we consider together the second 
and third terms of the l.h.s. of ([!]), defining the following 
Integral Structure Functions (ISF): 

D p (z, r) = ^ (Sv(z, r) 3 + ar ■ S{z) ■ Sv(z, r) 2 ) p/3 ) , (2) 

where a is an empirical prefactor of order one. These ISF 
are expected to take into account shear effects, particu- 
larly at large scales (see next paragraph) and to display a 
universal behavior. We consider here generic situations 
in which the shear reduces to S(z) = d z v x . Such situa- 
tion occurs near a rigid wall, where the principal mean- 
velocity component aligns in the x-direction, parallel to 
the boundary pjj| . The shear S(z) characterizes the vari- 
ation of v x along the z-direction, i.e. as one moves off the 
wall. Finally, we consider increments in the direction of 
the mean flow, i.e. orthogonal to the shear direction. 

ISF reduce on two different SF when either the first 
or the second term dominates. These two terms will ex- 
actly balance at scale L$ such that 5v{Ls)/Ls = a S. 
At scales r <C Ls, shear effects become negligible and 
homogeneous and isotropic scalings are expected. Kol- 
mogorov's scaling then yields Sv(r) ~ g 1 / 3 ^ 1 / 3 ^ where e 
denotes the mean energy dissipation rate. By extending 
this similarity relation to the scale r = L5, one obtains 
the usual dimensional estimate for the shear length scale 
Ls ~ |L2]]. In the logarithmic layer of a plane 

near- wall flow it yields L$(z) ~ z ||. Roughly speak- 
ing, Lg can be viewed as the size of small-scale eddies, 
whose turn-over time equals the shear time scale 1/5, 
imposed by the flow geometry and the stirring process 
of turbulence at large scales. Note that this estimate of 
the shear length-scale stems out from dimensional anal- 
ysis. In practice, there may be a prefactor in the expres- 
sion of L s : this is taken into account by the coefficint a. 
From previous reasoning it follows that D p (r) ~ D p (r) 
for r < L s and D p (r) ~ {rS) p/3 D 2p/3 {r) for r > L s . 

Another way to see that the central objects, in presence 
of shear, are the ISF, comes from the generalization of 
Yaglom's equation to homogeneous-shear Sows, i.e. with 



-er = D 3 (r) - %v h 

5 ar 

2s r 



dx x v' x (xq) ■ v' z (xq + x). 



(3) 



Supposing that this relation can be generalized, along the 
same line of idea of the Kolmogorov's Refined Similarity 
Hypothesis pM, one obtains 



Sv (r) + ar ■ S ■ v' x (x) • v' z (x + r) ~ e(r) • r, (4) 

where e(r) denotes the coarse-graining of the energy dis- 
sipation field e = f J2i (^iVj) , at scale r. 
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FIG. 2. In the turbulent boundary layer at the distance 
z + — 102. D3 and Dg (A) are compared with D3 and De (*). 
The scale r has been renormalized by the characteristic shear 
length-scale L s . The solid lines passing through D3 and D§ 
indicate the expected homogeneous and isotropic power-laws, 
respectively £3 = 1 and £g — 1.78. For comparison, the dashed 
line has slope 2. 

Simplifying, on pure dimensional grounds, the velocity- 
cross correlation v' x v' z , with Sv(r) 2 one ends up again with 
D p (r). Furthermore, we propose the following Refined 
Similarity Hypothesis 



D p {r) = /(6v(r) 3 + ar ■ S ■ 6v(r) 2 ) p/3 \ 



-( r)P /3\ r p/3_ 



This formulation is consistent, in the limiting cases of 
strong and negligible shear, with some recent findings 
(see ^, ||). In addition to that, the ISF should be able 
to abridge smoothly between these two limiting regimes, 
i.e. 

r«L 5 : D p (r) ^ (^e(r)P^y rP/ 3 (5) 
and r>i 5 : D p {r) ~ (e(r) p ^ . (6) 

Eqn. (||) is in agreement with the restoration of ho- 
mogeneity and isotropy at small scales. For r 3> Ls, 
eqn. (||) gives D 2 {r) ~ const (modulo possible loga- 
rithmic corrections), yielding for the energy spectrum 
E(k) ~ fc _1 (a relation suggested long time ago in Jl2|). 

In the previous picture, it is assumed that the shear 
length scale Ls remains larger than the dissipation length 
scale r\. The dissipation field e(r) is then expected to dis- 
play the same scaling properties as in homogeneous and 
isotropic turbulence. However, in regions where Ls iS r\ 
(very close to the wall) it is expected that shear effects 
act down to the dissipative scale and therefore can mod- 
ify the scaling behavior of e(r). The scalings of D p (r) 
should then also change according to (0) . 
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FIG. 3. From the boundary layer experiment, third-order 
structure function at various distances from the wall: z + = 37 
(o), z + = 124 (A) and z+ = 233 (□). The scaling properties 
of D3 (SF) do depend on the distance z + . On the contrary, 
Z?3 (ISF) displays the same scaling behavior for all z + . The 
dashed line has slope 1. The curves have been shifted verti- 
cally for convenience. 

The relevance of Integral Structure Functions for de- 
scribing the scaling properties of non-homogeneous shear 
flows is now tested on both experimental and numer- 
ical data. The experiment, performed in the recircu- 
lating wind tunnel of ENS-Lyon, consists in a turbu- 
lent boundary-layer flow over a smooth horizontal plate 
(see for details about the experimental apparatus). 
Velocity measurements are carried out at various eleva- 
tions from the plate in the logarithmic turbulent sub- 
layer Numerical results are obtained from a direct 
simulation of the Navier-Stokes equations in a rectangu- 
lar channel flow (see [Q for details). 

In Fig. ||, D3 and Dq, measured in the logarithmic 
boundary sublayer, are compared with the correspond- 
ing D3 and Dq. The shear S(z) = d z v x has been es- 
timated from the mean velocity profile. In standard 
non-dimensional variables namely z + — v^zjv and 
v + = v x /v* where v* is the characteristic velocity of the 
viscous sublayer, our data are well fitted by the loga- 
rithmic profile v + (z + ) = (1/k) log(z + ) + B. We obtain 
k ~ 0.4 and B ~ 5.26 in agreement with previously re- 
ported results |Q. For what concerns the coefficient a, 
all our results have been obtained with the fixed value 
a 0.2. Note that the constant a is not (a priori) in- 
tended to be universal but may strongly depend on the 
geometry and stirring process of the flow. However, it is 
expected to remain of order unity. In practice, the value 
of a has been extracted from data by requiring that one 
third-order ISF scale as r (see next paragraph). Corre- 
sponding estimates of L s are indicated in all figures. Fi- 
nally, one must point out that velocity increments have 
been estimated in the direction of the mean flow by use 
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FIG. 4. In the experimental boundary layer at z + — 102. 
SF and ISF, compensated by homogeneous and isotropic scal- 
ings, are displayed for p=l ... 6. 

of the Taylor hypothesis. Both D3 and Dq exhibit a 
power-law dependence on r but the scaling exponents 
are clearly different from those observed in homogeneous 
and isotropic (h-i) turbulence, respectively, C3 — 1 an d 
^6 ~ 1.78. On the other hand, the corresponding ISF 
exhibit power-laws in good agreement with h-i scalings 
(up to very large scales). 

In Fig. ||, third-order SF and ISF are displayed for 
various distances from the wall. We notice that scaling 
behavior of D3 changes with z + . On the contrary, all 
the corresponding D3 display the same power-law scal- 
ing with exponent 1. We recall that the coefficient a is 
kept constant and the shear is estimated from the mean 
velocity profile; there is no adjustable parameter. 

We now report a sharper test: SF and ISF, compen- 
sated by h-i power-law scalings, are displayed in Fig. |] 
for p = 1, . . . , 6 at distance z + = 102. A departure from 
h-i scalings is clearly observed for SF. 

On the contrary, D p (r)/r^ 1 ' exhibits a plateau up to 
very large scales, indicating that ISF roughly behave as 
r^ p . In other terms, ISF compensate shear effects and 
restore, via the extra term ar ■ S ■ Sv(r) 2 , the h-i scalings. 
Finally, the same test, made on numerical data, is re- 
ported in Fig. H for the sake of comparison. Data are ob- 
tained at distance z + — 25 from the wall, i.e. where the 
shear is strong. Results are in reasonable agreement with 
those of Fig. ||, despite the lower resolution of the numer- 
ical simulation. We would like to underline the major 
points of our study. Our description relates the scaling 
properties of velocity fluctuations in non-homogeneous 
shear flows, to those of the coarse-grained dissipation rate 
s(r). Provided that the shear length-scale remains much 
larger than the dissipative length-scale, we are inclined 
to believe that the scaling properties of e(r) remain flow- 
independent; the dissipation process, operating on very 
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FIG. 5. As in Fig. W but from numerical data at 



25. 



small-scales, is mainly insensitive to the presence of the 
shear. We then claim that the scalings of velocity fluctu- 
ations are different in sheared regions only because the 
similarity between Sv(r) and e(r) changes form (see (|^) 
and (^)). When L$ ^ T), shear effects acting down to 
dissipative scales are expected to modify the scaling be- 
havior of s(r). Nonetheless, we observe in Fig. [I] that 
relative scalings exponents of velocity structure functions 
remain universal. 

The introduction of ISF relies on quite simple dimen- 
sional arguments. However, they have proved to be valu- 
able (preliminary) tools in order to capture (at leading 
order) shear effects on the scaling behavior of velocity 
structure functions. We believe that some refinements 
in the definition of ISF may be possible, e.g. consider- 
ing the cross correlation v' x v' z instead of 5v(r) 2 : however 
ISF are of practical interest as they are easily accessible 
experimentally. 

A possible extension of ISF applies to thermal turbu- 
lent convection, the buoyancy term agST(r) playing a 
similar role of S ■ 6v(r). In that case, the ISF would read 



/3* B (r) = ((Sv(r) 3 + ar ■ agST(r) ■ Sv(r)) p/3 



and are expected to take into account buoyancy effects 
at scales larger than the Bolgiano length-scale. 

An important application of our findings concerns 
Large Eddy Simulations (LES). Usual eddy viscosity 
models are known to fail close to the boundary. As 
one moves near the boundary, the shear length-scale L$ 
becomes smaller and smaller. While L$ remains larger 
than the cut-off scale A of the LES, usual closure con- 
ditions, based on homogeneous and isotropic turbulent 
dynamics, remain acceptable. However, when L$ be- 
comes comparable or smaller than A, shear effects must 
be taken into account and the closure condition should 
be modified. An alternative consists in decreasing the 
mesh-size near the wall so that L$ always remains larger 
than A. More simply, our study suggests to consider 
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(5v(r) 3 + arSSv(r) 2 ) instead of 8v(r) in closure rela- 
tions. Along this line of idea, the Smagorinski's closure 
condition [l6| can be generalized in order to uniformly 
take into account shear effects. The eddy viscosity t'eddy 
then reads 

^eddy = C s A 2 ■ (S + aS) , 

where the extra term aS takes care of shear effects. C s is 
the empirical constant of the Smagorinski's closure and 
S denotes the rate of strain on scale A. 
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